Optimal multicopy asymmetric Gaussian cloning of coherent states 
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We investigate the asymmetric Gaussian cloning of coherent states which produces M copies from 
N input replicas, such that the fidelity of all copies may be different. We show that the optimal 
asymmetric Gaussian cloning can be performed with a single phase-insensitive amplifier and an 
array of beam splitters. We obtain a simple analytical expression characterizing the set of optimal 
asymmetric Gaussian cloning machines. 
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I. INTRODUCTION 

The perfect copying of unknown quantum states is for- 
bidden by the linearity of quantum mechanics Q . This 
observation lies at the heart of novel quantum commu- 
nication protocols, such as quantum key distribution 
(QKD) which allows the provably secure sharing of a 
secret key between two distant partners (see, e.g., Q). 
Any eavesdropping on a QKD system introduces noise 
into the transmission, which can be detected by the le- 
gitimate users. The optimal individual eavesdropping at- 
tacks on many QKD protocols consist of the optimal (ap- 
proximate) copying of the quantum states transmitted in 
the channel, where one copy is sent to the legitimate re- 
ceiver while the other copy is kept by the eavesdropper 
and measured upon at a later stage. After the seminal 
paper by Buzek and Hillery where the concept of 
universal quantum cloning machine was introduced, the 
issue of quantum copying has attracted considerable at- 
tention (see, e.g., [HOI). This effort culminated in the 
recent years with the experimental demonstration of op- 
timal 1 — > 2 cloning machines for polarization states of 
photons based either on parametric amplification |ES0 
or on the symmctrization of the multiphoton state on an 
array of beam splitters 0, 0] . The latter technique was 
also exploited lately to realize the universal symmetric 
cloning machine for qubits that produces three clones 

m 

In applications such as QKD, one is often interested in 
asymmetric cloning where the fidelities of the clones are 
different. This is indeed necessary to study the trade-off 
between the information gained by the eavesdropper and 
the noise detected by the legitimate users. The optimal 
1 — > 2 asymmetric cloni ng o f qubits and qudits has been 
studied in detail 0, IT3I [l4| and very recently an exper- 
imental demonstration of 1 — > 2 asymmetric cloning of 
polarization states of photons based on partial tele- 
portation |l6j | was reported. Going beyond two copies, 
multipartite asymmetric cloning machines have been in- 
troduced in |17|. which produce M copies with different 
fidelities Fj (j = 1, • • • M). Several examples of such mul- 
tipartite asymmetric doners for qubits and qudits were 
presented in 0, . 

In the context of the rapid development of quantum in- 



formation processing with continuous variables |2[J, the 
cloning of coherent states has been extensively studied 
over the last years 0, • It was shown that the opti- 
mal N — > M symmetric cloning of optical coherent states 
that preserves the Gaussian shape of the Wigner function 
can be accomplished with the help of a phase-insensitive 
amplifier followed by an array of beam splitters that dis- 
tributes the amplified signal into M modes [2^, 0] • It 
is also possible to exploit the off-resonant interaction 
of light beams with atomic ensembles and perform the 
cloning of coherent states into an atomic memory |25j . 
The cloning of a finite distribution of coherent states was 
studied |26J , and a reversal of cloning by means of local 
operations and classical communication was suggested in 



On the experimental side, the 1 — > 2 optimal Gaussian 
cloning of coherent states of light was recently demon- 
strated in [28|. There, the phase- insensitive optical am- 
plifier was replaced with a clever combination of beam 
splitters, homodyne detection, and feedforward, which 
effectively simulated the amplification process. Using ho- 
modyne detectors with very low electronic noise, it was 
possible to achieve a cloning fidelity of about 65%, very 
close to the theoretical maximum of 2/3. In another ex- 
periment, the 1 — > 2 telecloning of coherent states of light 
was also realized |29j . 

In this paper, we extend the concept of multipartite 
asymmetric cloning to continuous variables and present 
the optimal multipartite asymmetric Gaussian cloning 
machines for coherent states. These devices produce M 
approximate replicas of the coherent state |cv) from ./V 
input replicas, such that the fidelity Fj of each clone is 
generally different, and, for a given set of F\, . . . , Fm-i, 
the fidelity of the M-th clone Fm is the maximum possi- 
ble. The multicopy asymmetric 1 — * M cloning of coher- 
ent state was previously studied by Ferraro and Paris in 
the context of telecloning [3(j. Here, we rigorously prove 
that their scheme is optimal, and present a generic op- 
timal asymmetric cloning machine for any number N of 
input replicas, as well as its optical implementation. We 
also consider the related problem of the optimal partial 
state estimation of coherent states. 

The rest of the paper is structured as follows. In Sec- 
tion II, we present an optical cloning scheme based on 
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phase-insensitive amplification and passive linear optics. 
We also derive the trade-off between the fidelities (or, 
equivalently, added thermal noises), which fully charac- 
terizes the class of the optimal multipartite asymmetric 
Gaussian doners. In Section III, we describe an alterna- 
tive cloning scheme where the amplification is replaced by 
measurement and feedforward. In Section IV, we shortly 
discuss the relationship between the optimal asymmet- 
ric cloning and optimal partial measurement of coherent 
states. Then, the proof of the optimality of the asym- 
metric cloning machine is given in Section V. Finally, 
Section VI contains a brief summary and conclusions. 
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II. ASYMMETRIC GAUSSIAN CLONING OF 
COHERENT STATES 



FIG. 1: Optimal Gaussian N — > M fully asymmetric cloning 
of coherent states. See text for details. 



In what follows, we restrict ourselves to Gaussian 
cloning transformations. Note that it has been found 
in that the optimal 1 — > 2 cloning transformation for 
coherent states (i.e. the transformation that maximizes 
the single-clone fidelity) is in fact non-Gaussian, though 
the gain in fidelity is tiny. Nevertheless, it should be 
stressed that, while this non-Gaussian doner is optimal 
in terms of fidelity, it adds more noise to the clones (as 
measured by the quadrature variances) than the optimal 
Gaussian doner. In potential applications of cloning such 
as eavesdropping on QKD protocols with coherent states 
and homodyne detection [3^ H3 , one is often interested 
in minimizing the quadrature variance. In such a case, 
the Gaussian cloning turns out to be the most danger- 
ous attack (34[. Finding the multipartite generalization 
of asymmetric Gaussian cloning is therefore a very inter- 
esting question. 

As we will prove in Section V, the optimal Gaussian 
cloning machine has the simple structure depicted in 
Fig. 1, which is a direct generalization of the 1 — > 2 
asymmetric doner |24j . The signal contained in TV input 
replicas of the coherent state |a) is first collected into a 
single mode by an array of N — 1 unbalanced beam split- 
ters After this, a single mode a carries all the 
signal, and is in a coherent state \y/~Na). This mode is 
sent on an unbalanced beam splitter BS with amplitude 
transmittance t and reflectance r, which divides the sig- 
nal into two modes a and b\. Mode a is then amplified 
in a phase-insensitive amplifier (NOPA) with amplitude 
gain g. The modes a and b\ together with M — 2 aux- 
iliary modes bj, with j = 2, . . . , M — 1, are combined in 
a passive linear M-port interferometer IF whose output 
modes contain the M clones. The interferometer is de- 
signed in such a way that that the coherent component 
in each output mode is equal to a. In the Heisenberg pic- 
ture, the overall input-output transformation describing 
the doner depicted in Fig. 1 reads 
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where d is the creation operator of the idler port of the 
amplifier and bk are the annihilation operators of the 
M — 1 auxiliary modes, initially in the vacuum state. 
Here, rij represents the amount of noise added to the j- 
th clone, and the Kjk coefficients are chosen in such a way 
that the canonical commutation relations are conserved. 

It follows from the canonical transformations QJ that 
each clone is in a mixed Gaussian state with coherent 
amplitude a and added thermal noise characterized by 
a mean number of thermal photons rij. The Husimi Q- 
function of the j-th clone reads, 
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The fidelity of the j-th clone is proportional to the value 
of the Husimi Q-function at (3 = a, and is therefore a 
monotonic function of the added thermal noise, 
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The cloning is covariant and isotropic, i.e., the fidelity 
does not depend on the input state \a) and the added 
noise is the same for each quadrature. These are natu- 
ral conditions that the optimal cloning machine should 
satisfy. 

The shot-noise limited amplification is governed by the 
transformation 



= g(ta-rb 1 ) + y/g 2 



(4) 



and the total mean number of thermal photons produced 
during the amplification is n to t = 9 — 1- Since the linear 
interferometer does not add any noise, we have 



where nto 
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The total intensity of the coher- 
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ent signal after amplification is N(r + g t 



which 
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should be equal to Af|a| 2 if we require the coherent com- 
ponent of each clone to be equal to a. From this, we can 
determine the transmittance of BS, namely 



t = 



M-N 
"•tot N ' 



(6) 



The multiplets of rij cannot be arbitrary. Indeed, the M- 
port interferometer IF in Fig. 1 is described by a unitary 
matrix V, such that aj = YljLi 1 v jk°k + ^mow- The 
unitarity of V imposes a constraint on rij which can be 
expressed as 
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This formula provides a simple analytical parametriza- 
tion of the set of optimal N — * M multipartite asymmet- 
ric Gaussian cloning machines for coherent states. 

In the special case of a 1 — ► 2 asymmetric Gaussian 
doner, Eq. Q reduces to 



nmi = (1/2) 5 



(8) 



which coincides with the no-cloning uncertainty relation 
that was displayed in j23, 01 • (Note that 1/2 corre- 
sponds here to one shot-noise unit.) Interestingly, if wc 
consider a 1 — > 3 doner and assign to the first two clones 
the fidelity of the optimal 1 — *■ 2 symmetric doner, that 
is, n\ = ri2 = 1/2, we obtain by solving Eq. Q that the 
noise of the third clone is not infinite, = 2. As no- 
ticed in |18| , this means that some quantum information 
remains available beyond the one contained in the two 
clones (it actually corresponds to the information hidden 
in the anticlone). In the case where N = 1 but M is 
arbitrary, we recover the expression that was derived by 
Ferraro and Paris [3(j • Finally, note that if one clone is 
perfect, e.g., nu = 0, then Eq. J7J) is transformed into 
the same equation but for a (N — 1) — > (M — 1) doner, 
which means that one of the input replicas is simply redi- 
rected to the perfect clone while the cloning of the N — 1 
remaining input replicas into the M — 1 other clones is 
simply governed by the same relation. 



III. OPTIMAL ASYMMETRIC CLONING VIA 
MEASUREMENT AND FEEDFORWARD 



In the experimental demonstration of the optimal 
cloning of coherent states of light carried out in [28| . 
the amplification in a phase-insensitive amplifier was re- 
placed by a clever combination of a partial measurement 
and feedforward that effectively simulates the amplifica- 
tion process. Since an amplifier of arbitrary gain g can be 
implemented in this way [35| , the scheme shown in Fig. 1 
can be straightforwardly transformed into a setup which 
involves only passive linear optics, balanced homodyne 



FIG. 2: Setup for the multipartite asymmetric Gaussian 
cloning of coherent states using homodyne detection and feed- 
forward. 



detection and coherent displacement of the beams pro- 
portional to the measurement outcomes. The resulting 
configuration is illustrated in Fig. 2. 

We assume that all available signal has been collected 
into a single mode a which is thus in the coherent state 
\yNa). The beam is divided into two parts on a beam 
splitter BS g with amplitude transmittance t and re- 
flectance f. The reflected part is fed into a heterodyne de- 
tector consisting of a balanced beam splitter whose auxil- 
iary input port c is in the vacuum state and two balanced 
homodyne detectors BHD measuring the x and p quadra- 
tures, respectively. This detector effectively measures the 
operator o = r a + tb+c^ . The portion of the beam trans- 
mitted through BS 5 , characterized by la — fb, is divided 
into M modes aj by an array of M — 1 beam splitters 
with reflectances r k . Each mode a k is then coherently 
displaced by amount gj o, where gj is the electronic gain 
of the corresponding feed forward. The added thermal 
noise in the output mode aj reads rij = glj, which imme- 
diately fixes all electronic gains, 



9j 



(9) 



As shown in |35|. the optical amplification gain g is ob- 
tained in the feed-forward equivalent scheme if the beam 
is split on a beam splitter with reflectance yl — 1 /g 2 , 
and the reflected part is heterodyne measured. Since 
in the scheme of Fig. 1 only an (amplitude) fraction t 
of the input beam a is actually amplified, we see that 
by replacing the amplifier with a beam splitter of re- 
flectance y/l — 1/.9 2 , the fraction of the beam that is 
sent to the heterodyne detector is t yl — l/<? 2 - Thus, 
in Fig. 2, the reflectance of the beam splitter BS 9 must 
be r — t y 7 ! — 1/g 2 , which yields 
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It remains to determine the reflectances rk of the final 
array of beam splitters. They are fixed by the condition 
that the coherent amplitude of each clone is a. After 
some algebra we find that 



VI + ntot - y/(M - N) rij _ ^-1/2 



^{2 + n tot )N-M 



(11) 
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From this formula, all r^'s can be calculated in an itera- 
tive way, starting from r*i [which is given by Eq. for 
j = 1 with the product over k replaced by 1]. 



IV. OPTIMAL CLONING AND OPTIMAL 
PARTIAL ESTIMATION OF COHERENT STATES 

There is a close relationship between optimal cloning 
and optimal state estimation. An interesting scenario 
that recently attracted a lot of attention consists in the 
partial estimation of a state, which yields the classical 
estimate of the state as well as the perturbed quantum 
state S El Ell El El According to the fact that 
in the limit of an infinite number of copies, the optimal 
cloning becomes equivalent to optimal state estimation 
|36j |. this optimal partial estimation can be viewed as 
a limiting case of an asymmetric cloning producing one 
(quantum) copy with fidelity F and infinitely many (clas- 
sical) copies with fidelity G ^t|- From the analytical 
formula (JJJ, we can thus rigorously derive the optimal 
trade-off between the fidelities F and G in the partial 
Gaussian estimation of coherent states. We set ni = rip, 
Hi = Tis = . . . = njvf = tig and take the limit M — > 00, 
which results in 



(n F + l) 2 



(12) 



In the limit of an undisturbed quantum copy [rip = 0), 
we have an infinitely noisy state estimation, as expected. 
We also note that np = uq = 1 is a solution of Eq. (|T2*)l . 
which corresponds to the optimal (full) estimation of co- 
herent states. Equation (|12[) also translates in the fol- 
lowing relation between the fidelities, 



G = 



4F(1 - F) 
4F(1-F) + 1' 



(13) 



This agrees with the trade-off derived in (33 which con- 
firms that the experimentally demonstrated partial mea- 
surement of coherent states in that work was indeed op- 
timal among all Gaussian strategies. Note that using a 
non-Gaussian protocol a slightly better trade-off between 
F and G could be achieved 



fitly 



V. PROOF OF OPTIMALITY 

In what follows, we will prove the optimality of the 
asymmetric doner defined in Sees. II and III. Let us first 



note that the fidelities are monotonic functions of rij, 
so that instead of maximizing the fidelities Fj we can 
equivalently minimize the added thermal noise rij. The 
design of the optimal doner can be thus rephrased as the 
minimization of a cost function p"?| 



M 



(14) 



which is a linear convex mixture of rij, Xj > 0. The 
ratios of the coefficients Xj control the asymmetry of the 
cloning machine. 

The most general Gaussian operation is a trace- 
preserving Gaussian completely positive (CP) map [2l| . 
and we have to minimize l|14|) over all such maps. At the 
level of covariance matrices 7, the Gaussian CP map acts 
as 



7out = Sji n S + G. 

The covariance matrix of N modes is defined as 



(15) 



(ArjArk + Ar^Arj), where 



Ijk = 
, xn,Pi ■ ■ ■ ,Pn) 
is the vector of quadrature operators, [2j,pfc] = idjk- The 
first moments transform under the Gaussian CP map ac- 
cording to (r ou t) = S(r in ). 

The matrices S and G must satisfy the complete posi- 
tivity constraint 



A = G + iK > . 



where the matrix 



K = J Mo 



I 

-I 



S J M in S 



(16) 



(17) 



comprises the commutators of the quadrature operators, 
while / denotes the identity matrix of dimension M, and 
Mi n and M out are the number of input and output modes, 
respectively. The cloning machine of interest has effec- 
tively only a single input mode a (as we collect the N 
input signals into a single mode) and M output modes, 
so that M; n = 1 and M out = M. The M x 2 matrix S is 
fixed by the condition that the first moments should be 
preserved by cloning. We get 



S T = 



1 



1 1 





1 
1 1 



(18) 



We have to minimize C(rij) over the set of all Gaussian 
completely positive maps ()15|l with the matrix S given by 
l)18|l. that is, we have to optimize over all G's satisfying 
(|16|l . Clearly, if G\ and G2 satisfy (|16|) . then any convex 
combination pG\ + (1 — p)G2 with p S [0, 1] also does. 
Thus, we have a convex optimization problem. Moreover, 
since rij = (Gjj + GM+j,M+j + 2/N — 2)/4, the cost 
function (|14f) is linear in the matrix elements of G. Hence, 
the problem amounts to minimizing 



C(G) 



M 

E 



x 3\Gii + GM+j,M+j) 



(19) 
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under the constraints (|16H , which is an instance of a linear 
scmidefinite program |4J|. We shall now prove the opti- 
mality of JQ) by deriving a lower bound on G(G) which 
is saturated by (JTJ. 

The specific feature of the transformation is that 
only a single creation operator cr is admixed to the an- 
nihilation operators. This operator is responsible for the 
added noise in cloning. Since all modes are initially in 
coherent states, the normally ordered moments of the 
operators a,j for the clones can be easily calculated, 



(A aj Aa k ) = (AatA<4) = 
(AajAafc) = y/njnk. 



0, 



The covariance matrix of the M clones is then fully de- 
termined by the added noises rij, 



Tout 



I + 2F 
I + 2F 



(20) 



where F is a symmetric M x M matrix with elements 
Fjk = y/rijUk- Since the matrix S is fixed and the input 
state of the doner is a coherent state with covariance 
matrix 7; n = I, the matrix G op t corresponding to trans- 
formation can be determined from Eq. i|15|) . This 
yields 



G 



opt 



2F-±H 






I+2F-±H 



(21) 



where H is a matrix whose elements are all equal to one, 
Hjk = 1. 

Since the transformation Q can be associated with 
a CP map, the matrix G op t must satisfy the inequality 
pfijl. For the particular S 1 matrix (|18|1 . this gives 



_/7 + 2F-iff i(I-iff) 

^P* - I .Vr 1 m r i oc 1 



-i(I-±H) I + 2F 



N 



H 



> (22) 



We can transform the matrix A opt to a block-diagonal 
form with the help of the unitary matrix 



U 



( 



I il 



y/2 \U I 



(23) 



which gives UA opt W = diag(2,F + 21 - 2H/N,2F), so 
that there remains to prove that F + I — H/N > and 
F > 0. The matrices F and H both have rank one, 
and can be written in Dirac notation as F = |/)(/|, 

where |/) = EjliV™^')' and H = whcrc 
\h) = Ylj=i li)- The condition F > is satisfied by defi- 
nition, while solving F+I—H/N > yields the nontrivial 
constraint on n^s. Actually, it is sufficient to check the 
latter positivity condition in the two-dimensional sub- 
space spanned by the (un-normalized) vectors |/) and 
| ft), and one can prove that it is indeed satisfied provided 
that Eq. holds. 



We now derive a tight lower bound on C(G), which is 
saturated by (JIJ . Suppose that we find a positive semidef- 
inite matrix Z > such that it satisfies the conditions 



and 



Tr[ZG] = G(G) 



ZA opt = Z{G opt + iK) = 0, 



(24) 



(25) 



where iK — iJM — iSJ%S T . Then, the Gaussian CP map 
with matrix G op t is the optimal one that minimizes G (G) . 
Since for every admissible G we have G + iK > 0, it 
follows from Z > that Tr[Z(G+iK)] > 0, which implies 
that G(G) > -iTi[ZK], VG. Equation (gSJ implies that 
this lower bound is saturated by G op t , which is therefore 
optimal. 

The matrix Z can be determined from Eqs. (|24|) and 
G3. We find that 



Z = 



X iY 
-iY X 



(26) 



where X = diag(xi, . . . ,xm) is fixed by Eq. ()24fl. while 
Y is a real symmetric matrix that satisfies 

Y(I - N^H) + X(I + 2F- N~ l H) = , 
XF — YF , (27) 

as a consequence of Eq. I|25|l . Note that X > by def- 
inition because Xj > 0. Since the matrix I — N^ 1 !! is 
invertible, we can express the matrix Y in terms of X 
using the first condition of (|2*7|) . 



Y = -X[I + 2F(I - (M - N)- 1 !!)] 



(28) 



The second condition of l|27|) is then satisfied for any X 
provided that J7J) holds. 

In order to further simplify the matrix Y, we need to 
establish the relationship between Xj and nj. Without 
loss of generality, we can restrict ourselves to the cloning 
machines that satisfy l[7)l. and minimize the cost C(rij) 
under the constraint Q. Using the standard method of 
Lagrange multipliers, we obtain the extremal equations 
for the optimal rij's for a given set of Xj's, 



A I 



l JV 3 



X(M — N)y/nj + y^nk = 0. 



(29) 



k=l 



with A being the Lagrange multiplier. With the help of 
these formulas, we can show that 



XF(I - (M -N)- 1 !!) 



1 



X(M - N) 



XFX, 



so that 



Y = -X 



X(M - N) 



XFX. 



(30) 



(31) 
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The matrix XFX is symmetric, hence Y = Y T as re- 
quired. 

The last step of the proof is to show that the matrix Z 
is positive semidefinite. We first apply a transformation 
that preserves the positive semidefmiteness, Z = VZV', 
where V = diagpf" 1 / 2 ,X~ 1 / 2 ), 

-_( I -il -i2r l X 1 / 2 FX 1 / 2 \ 

Z ~ \ il + i2 V X 1 / 2 FX 1 / 2 I J 1 

where 77 = 1/(A(M - N)). We multiply Eq. (J2U with 
X~ x and take the trace, so we find that 

v Tt[X 1/2 FX 1/2 } = -1, (32) 

where we made use of Eq. (JJJ. Since F is propor- 
tional to rank one projector the normalization (|32|) im- 
plies that rjX^-^FX 1 / 2 = -\<t>){(f>\ = where \<f>) 
is a normalized real vector, (4>\4>) = 1. We can con- 
vert Z to block diagonal form with the unitary 1)23(1 . 
UZW = diag(2<I>, 21 - 2$), which is obviously positive 
semidefinite. This concludes the proof of optimality of 
the cloning machine Q). 

VI. CONCLUSIONS 

In summary, we have proposed a multipartite asym- 
metric Gaussian cloning machine for coherent states. The 



machine produces M approximate copies of a coherent 
state from N replicas of this state in such a way that 
each copy can have a different fidelity. A simple ana- 
lytical formula characterizing the set of optimal Gaus- 
sian asymmetric cloning machines has been derived, and 
it was shown that the asymmetric cloning can be real- 
ized by amplifying of a part of the input signal followed 
by mixing the amplified signal and the bypass signal to- 
gether with auxiliary vacuum modes on an array of beam 
splitters with carefully chosen transmittances. An alter- 
nate implementation is also described, where the ampli- 
fier is replaced by a passive optical circuit supplemented 
with feedforward. We hope that our study of multipartite 
asymmetric cloning will trigger further investigations of 
optimal quantum information distribution in continuous- 
variable quantum communication networks. 
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